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Abstract 

This note aims at demonstrating the advantage of moving-water well-balanced 
schemes over still-water well-balanced schemes for the shallow water equations. We 
concentrate on numerical examples with solutions near a moving-water equilibrium. 

For such examples, still-water well-balanced methods are not capable of capturing the 
small perturbations of the moving-water equilibrium and may generate significant spu¬ 
rious oscillations, unless an extremely refined mesh is used. On the other hand, moving- 
water well-balanced methods perform well in these tests. The numerical examples in 
this note clearly demonstrate the importance of utilizing moving-water well-balanced 
methods for solutions near a moving-water equilibrium. 
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1 Introduction 


The main objective of this note is to demonstrate the advantage of schemes which are well 
balanced for moving water equilibrium, over those which are only well balanced for still water, 
for small perturbations of the moving water equilibrium for the shallow water equations. 

Well-balanced schemes refer to those schemes which have zero truncation error for certain 
steady-state solutions of hyperbolic balance laws 

Ut + f{U), = g{U,x). 

That is, for certain non-trivial functions V(x) satisfying 

f{V), = g{V,x), 

the well balanced schemes do not have any truncation error. Note that in general V (x) is 
unknown and is not a polynomial, thus to require a zero truncation error is a difficult task. 

For shallow water equations, there are typically two classes of well-balanced schemes, 
namely those which are well balanced for still-water equilibrium, where the velocity is zero, 
and those which are well balanced for general moving-water equilibrium. A scheme which 
belongs to the latter class is more difficult to construct. We refer to [1] for a recently 
developed high-order accurate finite volume weighted essentially non-oscillatory (WENO) 
scheme which is well balanced for general moving-water equilibrium. Discussions on the 
history of the development of well balanced schemes for shallow water equations, with an 
extensive reference list, can be found in [2]. 

Since it is much more difficult to construct well-balanced schemes for moving-water equi¬ 
librium than for still water, a natural question is whether the former has any advantage over 
the latter. However, to our best knowledge, examples to show this advantage seem to be 
absent in the literature. In the Castro Urdiales meeting for which this special issue is dedi¬ 
cated to, there was an intense discussion on this issue but no such examples were provided 
in the discussion. 
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The purpose of this note is therefore to provide carefully selected numerical examples 
to demonstrate the advantage of well-balanced schemes for moving-water equilibrium over 
well-balanced schemes for still water. These examples have solutions which contain small 
perturbations from a moving-water equilibrium. A well-balanced scheme for still water would 
generate spurious oscillations at the level of the truncation error of the scheme, until the 
mesh is extremely rehned. On the other hand, a well-balanced scheme for the moving-water 
equilibrium could resolve these small perturbations on much coarser meshes. These examples 
clearly demonstrate the advantage of well-balanced schemes for moving-water equilibrium. 
We shall use the high order hnite volume WENO schemes in [1], which are well balanced 
for still water, and those in |T], which are well balanced for moving water, in the numerical 
tests. 

2 Shallow water equations and well-balanced methods 

In one space dimension, the shallow water equations take the form 



( 2 . 1 ) 


where h denotes the water height, u is the velocity of the fluid, b represents the bottom 


topography and g is the gravitational constant. The still-water (also referred as lake at rest) 
steady state is given by 


M = 0 and h + b = const, 


( 2 . 2 ) 


which represents a still flat water surface. The general moving-water steady state is given 
by 



(2.3) 


We introduce the notation U = {h, hu)'^ for the conservative variables, with the superscript 
T denoting the transpose, and V = {m, EY' for the equilibrium variables. We refer to [1] 
for the variable transformation between U = U{V,x) and V = V{U,x). 
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The computational domain is discretized into cells Jj = i = 1, ■ ■ ■ , N. We 

denote the size of the i-th cell by Axi and the center of the cell by x* = + ^*+1 j /2- 

Let U{xi,t) = Jj U{x,t) dx/Axi denote the cell average of U{-,t) over the cell Jj. The 
computational variables are f/”, which approximate the cell averages U{xi,t) at the time 
t = For the ease of presentation, we denote the shallow water equations fl2.ip by 

Ut + f{U)^ = g{h,b). 

2.1 Methods preserving still-water equilibrium 

Many high-order well-balanced methods for the still-water steady state solution fl2.2p have 
been developed in the literature for the shallow water equations, see for example [2] for a list 
of references. In this note, we will use the hfth-order hnite volume WENO methods based 
on the separation of the source terms developed in Xing and Shu [1]. In this subsection, we 
briefly review this well-balanced scheme in one dimension. We refer to |1] for further details. 

We hrst reconstruct the point values i at the cell interface from the given cell averages 
Ui by the standard WENO reconstruction procedure, which can be eventually written out 
as 

r r—1 

wuU^+k = S^{U)i, 17-1 = AktJ^+k = S^{U), (2.4) 

k=—r-\-l k=—r 

where r = 3 for the hfth order WENO approximation and the coefficients Wk and Wk depend 
nonlinearly on U. Now we freeze the cell averages U and the corresponding coefficients 
Wk and Wk- With these hxed coefficients, we apply the linear operators W i—>• S^{W) to 
W = (fe, 0)^ to compute the reconstructed values Assume U is the still-water steady 

state solution satisfying fl2.2p . we clearly have 

h^.+b^i= constant. (2.5) 

The main idea in constructing well-balanced methods is to decompose the integral of the 
source term into a sum of several terms, then compute each of them in a way consistent 
with the approximation for the corresponding flux terms. We hrst rewrite the shallow water 
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equations fl2.ip as 


ht + {hu)^ = 0 

{hu)t+ + -g{h + h)h^, 

for which the semi-discrete numerical method takes the form 


( 2 . 6 ) 


= -(F,+i - h-i) + G(u, x), 


(2.7) 


where F is a numerical flux, such as the Lax-Friedrichs flux 


F, 1 = 


*+ 


mzi) + f(Gii)-a 


h + h 

hu 


h + b 
hu 


i+2 V / i+ 


( 2 . 8 ) 


with a = max[/ \u ± \/^|. The first component of the source term approximation G{U,x) 
is zero, and the second component is given by 

(''i+1 “ ''i-i) -s('‘ + *)i(4+5 -Si-j) ~ (ft + * - (A + ^i) hdx, (2.9) 


w.th W,. = i ((ft.;p^ + (ft.:,aud 6„. = I(ft+. + ft.:, J. 


2.2 Methods preserving moving-water equilibrium 


We have developed high-order well-balanced hnite volume WENO schemes for the moving- 
water equilibrium (12.dh for the shallow water equations in [1]. In this subsection, we briefly 
review these methods and refer to [1] for further details. 

At each time step, we hrst apply the usual WENO reconstruction procedure to the cell 
averages f/j, and obtain hence 

Given the cell averages Ui and a bottom function b{x), we choose local reference values 
Vi of the equilibrium variables. These are defined implicitly by the requirement that 

^ f U{V,x)dx = Ui. (2.10) 

I J II 

Relation (12.10^ chooses Vi as the unique (see the paragraphs preceding [H Def.3.2]) local 
equilibrium such that the corresponding conserved variables U{Vi,b{x)) have the same cell 
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average Ui as the numerical data. It is proven in PQ Def.3.2] that, if the data U{x) and h{x) 

are in local equilibrium (y{U{x),x) = V for all cells Jj), the reference equilibrium states 

Vi computed via (I2.inp coincide with the true local steady state V. The reconstruction is 

completed by limiting the reconstruction with respect to the reference values Vi (see [H 

(3.18)]), to obtain the limited values and V. 

1 + 2 

The fourth-order well-balanced scheme is given by 


Ax, 


dU{xi, t) 
dt 


=-(F(£/-,, £/+.)-Fie/-,, £ 7 +, )) + s., 


( 2 . 11 ) 


where the function is a conservative, Lipschitz continuous numerical flux, and 







( 2 . 12 ) 


The total source term Si is given by 

AQ — Q 

:= + fiUlA - fiUlA + /(g^i) - /(C-J, (2.13) 

where fVi = f^(C^i,6^i). The extrapolated interior source term (45*2 — *S'i)/3 is dehned 
by 



and the well-balanced quadrature of the source term s*"’* is given by 

sT\Ul, Ur, bL, bn) = -^g{hL + hR){bR - b^ + sf* 


where s)"* is given by [H (3.63)-(3.65)]. 


(2.14) 

(2.15) 


(2.16) 


3 Numerical examples 

In this section we present numerical results of both moving-water well-balanced methods and 
still-water well-balanced methods presented in Section [2l for a carefully selected set of test 
examples in one dimension. Comparison of these results is provided as a demonstration of the 
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advantage of moving-water well-balanced methods. In all the examples, time discretization 
is by the classical third-order total variation diminishing (TVD) Runge-Kutta method [3], 
and the CFL number is taken as 0.6. The gravitation constant g is taken as 9.812m/s^. 


3.1 Perturbation of a moving water equilibrium 

The following test cases are chosen to demonstrate the capability of these schemes for com¬ 
putations on the perturbation of steady state solutions. 

The bottom topography is given by: 



0.2-0.05(a:-10)2 if 8 < x < 12, 
0 otherwise. 


(3.1) 


in the computational domain [0,25]. Three steady states, subcritical or transcritical flow 
with or without a steady shock will be investigated. 

Our initial conditions are given by imposing a small perturbation of size 0.05 on the height 
of these steady states in the interval [5.75, 6.25]. Theoretically, this disturbance should split 
into two waves, propagating to the left and right respectively. Note that in [1] , we have shown 
the numerical results of the moving-water well-balanced methods with a perturbation size 
0.01, which demonstrate that these small perturbations are well captured. Here we provide 
the numerical results of both the moving-water well-balanced methods and still-water well- 
balanced methods, and demonstrate the different behavior of these two methods, 
a): Subcritical flow: 

The initial condition is given by: 


(3.2) 


E = 22.06605, m = 4.42, 


together with the boundary condition that the discharge m=4.42 is imposed at upstream 
and the water height h=2 is imposed at downstream when the flow is subcritical. 

We run the tests with 100 uniform cells until t = 1.5, when the downstream-traveling 


water pulse has already passed the bump. The differences between the water height h at 


that time and the background moving water state, are shown in Figure 13.11 The result of 
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Figure 3.1: The difference between the height h at time t = 1.5 and the background moving 
steady state water height fl3.2p . when 100 uniform cells are employed. An initial perturbation 
of size 0.05 is imposed between [5.75, 6.25]. Left: result based on still-water well-balanced 
scheme. Right: result based on moving-water well-balanced scheme. 

the still-water well-balanced method is plotted on the left and that of the moving-water well- 
balanced method is shown on the right. The difference can be easily observed, as there are 
significantly more spurious oscillations generated by the still-water well-balanced method. 
As we refine the mesh to 1000 cells, the results are shown in Figure [321 Because Ax is now 
very small, the perturbation is relatively big in comparison with truncation errors of the 
schemes and can be well captured by both methods, 
b): Trans critical flow without a shock: 

The initial condition is given by: 

1 53 ^ 

E = -^- r + 9.812 X 0.66, m = 1.53, (3.3) 

2 X 0.66^ 

together with the boundary condition that the discharge m=1.53 is imposed at upstream 
and the water height h=0.66 is imposed at downstream when the flow is subcritical. 

We run the tests with 100 uniform cells until t = 1.5, when the downstream-traveling 
water pulse has already passed the bump. The differences between the water height h at that 
time and the background moving water state, are shown in Figure 1321 Again, the difference 
can be easily observed, as there are significantly more spurious oscillations generated by the 

















Figure 3.2: Same as in Figure 13711 but 1000 uniform cells are employed. 


still-water well-balanced method. As we refine the mesh to 1000 cells, the results in Figure 
13.41 show resolved solutions for both methods, 
c): Transcritical flow with a shock: 

The initial condition is given by: 


E = 


:(9.812 X 0.18)3) + 9.812 x 0.2 if x < 11.665504281554291 


0.18" 


+ 9.812 X 0.33 


otherwise 


m = 0.18, (3.4) 


2 X 0.33" 

together with the boundary condition that the discharge m=0.18 is imposed at upstream 
and the water height h=0.33 is imposed at downstream. 

We run the tests with 200 uniform cells until t = 3. The differences between the water 
height h at that time and the background moving water state, are shown in Figure 13.51 
Once again, the difference can be easily observed, as there are signihcantly more spurious 
structures generated by the still-water well-balanced method. As we refine the mesh to 1000 
cells, the results are shown in Figure 13^ Even at this resolution there is still an advantage 
using the moving-water well-balanced method. Note that near the point x = 11.7 there is a 
big displacement for both methods, as we recall that 11.7 is the position where the shock is 
located. 
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Figure 3.3: The difference between the height h at time t = 1.5 and the background moving 
steady state water height fl3.3p . when 100 uniform cells are employed. An initial perturbation 
of size 0.05 is imposed between [5.75, 6.25]. Left: result based on still-water well-balanced 
scheme. Right: result based on moving-water well-balanced scheme. 




Figure 3.4: Same as in Figure 1331 but 1000 uniform cells are employed. 
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Figure 3.5: The difference between the height h at time t = 3 and the background moving 
steady state water height fl3.4p . when 200 uniform cells are employed. An initial perturbation 
of size 0.05 is imposed between [5.75, 6.25]. Left: result based on still-water well-balanced 
scheme. Right: result based on moving-water well-balanced scheme. 




Figure 3.6: Same as in Figure 1331 but 1000 uniform cells are employed. 
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3.2 Perturbation with smaller magnitude 

In this subsection, we utilize a smaller perturbation to these tests and show that moving- 
water well-balanced methods demonstrate more clearly their advantage in capturing such 
perturbations. 

We keep the main setup the same as in Section 13.11 and impose a smaller perturbation of 
size 0.001 to these steady states. To save space, only the results of the subcritical flow (13.2p 
are shown. The differences between the water height h at that time and the background 
moving water state, when 100 uniform cells are used, are plotted in Figure ITTH We also show 
the differences of the momentum hu in Figure [T8l These figures clearly demonstrate that the 
still-water well-balanced method is not capable of capturing such a small perturbation on the 
coarse mesh, as large spurious oscillations are observed. As we rehne the meshes to 200 cells, 
the results are shown in Figures [3.91 and [3.101 where spurious oscillations for the still-water 
well-balanced method have a reduced magnitude, but its performance is still significantly 
inferior to that of the moving-water well-balanced method. The results with 1000 cells are 
shown in Figures [3.Ill and 13.121 Even though the difference between the two methods is now 
signihcantly reduced, the advantage of the moving-water well-balanced method can still be 
observed on such a rehned mesh. 

Acknowledgement. The hrst author is a contractor [UT-Battelle, manager of Oak 
Ridge National Laboratory] of the U.S. Government under Contract No. DE-AC05-00OR22725 
Accordingly, the U.S. Government retains a non-exclusive, royalty-free license to publish or 
reproduce the published form of this contribution, or allow others to do so, for U.S. Govern¬ 
ment purposes. 

References 

[1] S. Noelle, Y. Xing and C.-W. Shu, High order well-balanced finite volume WENO 
schemes for shallow water equation with moving water, J. Comput. Phys. 226 (2007), 


12 

















Figure 3.7: The difference between the height h at time t = 1.5 and the background moving 
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Figure 3.8: The difference between the momentum hv at time t = 1.5 and the background 
moving steady state water height fl3.2p . when 100 uniform cells are employed. An initial 
perturbation of size 0.001 is imposed between [5.75, 6.25]. Left: result based on still-water 
well-balanced scheme. Right: result based on moving-water well-balanced scheme. 
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Figure 3.9: Same as in Figure ETl but 200 uniform cells are employed. 




Figure 3.10: Same as in Figure |321 but 200 uniform cells are employed. 
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Figure 3.11: Same as in Figure ETl but 1000 uniform cells are employed. 




Figure 3.12: Same as in Figure |3S1 but 1000 uniform cells are employed. 
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